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ABSTRACT 

Aims. We present a study of the distribution of stellar ages in the Orion Nebula Cluster (ONC) based on accurate HST photometry 
taken from the HST Treasury Program observations of the ONC utilizing the most recent estimate of the cluster's distance (Menten et 
al. 2007). We investigate the presence of an intrinsic age spread in the region and a possible trend of age with the spatial distribution. 
Methods. We estimate the extinction and accretion luminosity towards each source by performing synthetic photometry on an em- 
pirical calibration of atmospheric models (Da Rio et al. 2010) using the package Chorizos (Maiz-Apellaniz 2004). The position of 
the sources in the HR-diagram is compared with different theoretical isochrones to estimate the mean cluster age and age dispersion. 
Through Monte Carlo simulations we quantify the amount of intrinsic age spread in the region, taking into account uncertainties on 
the distance, spectral type, extinction, unresolved binaries, accretion and photometric variability. 

Results. According to Siess et al. (2000) evolutionary models the mean age of the Cluster is 2.2 Myr with a scatter of few Myrs. With 
Monte Carlo simulations we find that the observed age spread is inconsistent with a coeval stellar population, but is in agreement with 
a star formation activity between 1.5 and 3.5 Myrs. We also observe light evidence for a trend of ages with spatial distribution. 

Key words. Stars: formation - Stars: low-mass - Stars: pre-main sequence - open clusters and associations: individual: Orion Nebula 
Cluster - Methods: statistical - Techniques: photometric 



1. Introduction 

The analysis of the age distribution of pre-Main Sequence (PMS) 
stars in young clusters represents the most direct way of deduc- 
ing how star formation proceed s in molecular clouds dHartmannl 
I l2001t iHillenbrand etall 120081: iDa Rio et alJ l2010bl) . Ages of 
young stellar objects are generally estimated from the compari- 
son between the position of individual stars on the HR-diagram 
and theoretical isochrones and evolutionary tracks. For the ma- 
jority of stars, the main uncertainties on age are related to the 
determination of the intrinsic luminosity, since the pre-main- 
sequence evolution of low mass stars takes p lace at nearly con - 
stant effective temperature. As pointed out by Hartmann (2001), 
the main sources of error on the luminosity are relative to the 
source distance, photometric variability, extinction, unresolved 
binary companions, and accretion. The discrepancies between 
different families of pre-main-sequence evolutionary models 
add further uncertainty, in particular the different slopes of the 
isochrone lines in the HR diagrams make it difficult to unequiv- 
ocally derive not only the absolute but also the relative ages of 
cluster members. 

Moreover, other physica l processes may a ffect the intrinsic 
luminosity of PMS stars. iBaraffe et alJ (l2009h present calcula- 
tions where episodic accretion at early stages of evolution of 
low-mass stars provides a possibl e explanation for t he observed 
spread in the HR diagram. iHosokawa et alJ d201 ll) show that, 
even though episodic accretion may occur, it is not able to ex- 
plain the observed scatter. In fact, the spread obtained varying 
the initial radius or the thermal efficiency is much larger than the 
spread that results from different accretion histories, even though 
it is not enough to account for all the observed spread in the HR 



diagram ( IHosokawa et alj201 lb . However, on one hand, a model 
quantifying whether the accretion-induced luminosity spread is 
able to account for the observations has not been created yet. On 
the other, we lack a direct observational proof for this scenario. 
Therefore we still interpret the distribution of radii for a given 
mass as a spread in age. 

In this paper we revisit this problem on the basis of new 
data from the Orion Nebula Cluster obtained with the Hubble 
Sp ace Telescope (HST). The same issue was partially discussed 
by IDa Rio et alJ d2010al) . using ground based photometry. Our 
present goal is to assess the possibility of inferring robust age 
distributions by using high precision multi-band photometry. 
Fitting all available photometric data at the same time, we can 
have accurate estimates for extinction and accretion luminosity 
which allows us to better constrain the errors on the luminosity 
and the age of each source. We also focus on the assessment of 
apparent age spread induced by these errors. After a brief de- 
scription of the photometry used and the selection criteria of our 
source sample (Sect. |2|, we present the analysis carried out to 
assemble an observed HR diagram for the ONC (Sect. [3). The 
study of the relative uncertainties and the results obtained are de- 
scribed in Sect. |4] and Sect. [5] Finally, we discuss the possibility 
of a true age spread in the Cluster (Sect. |6]l and the dependence 
of age on the spatial distribution (Sect. [7]i- 

2. Photometry 

The data used for the present work consist mainly of broad-band 
photometry extracted from the large data-set of ACS/WFC ob- 
servations collected for the HST Treasury Program on the Orion 
Nebula Cluster (GO-10246, PI M. Robberto). This program has 
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Table 1. ACS/WFC filters 



Filter 


Ground Equivalent 


Exposure time (s) 


F435W 


Johnson B 


420 


F555W 


Johnson V 


385 


F658N 


Ha +[N II] ,16583 


340 


F775W 


Cousin IC 


385 


F850LP 


z-Band 


385 



covered an area of about 450 arcmin 2 , roughly centered on the 
Trapezium Cluster in five filters. In Table [TJ we list the ACS fil- 
ters, their ground-based telescope equivalents and the duration 
of individual exposures. The complete ACS/WFC catalog for 
the ONC includes ~3200 sources (Robberto 2011, in prepara- 
tion) that have been classified as either a single star, binary star, 
photo-evaporated disk, dark silhouette, or possible background 
galaxy. A subset of th e catalog de s cribin g the non stellar objects 
has been presented bv lRicci et al.l (120081) . 

The deep ACS observations, designed to sample the substel- 
lar component of the cluster, do not provide photometry for the 
brightest sources because of saturation. We have therefore com- 
plemented the ACS survey with ancillary ground based data, 
collected within the same HST Treasury Program. These ob- 
servations were obtained simultaneously at the MPG-ESO 2.2- 
m Telescope at La Silla with the Wide Field Imager (WFI) 
dDa Rio et al.l l2009f) and with the Infra red Side Port Image r 
(ISPI) at the CTIO Blanco 4-m Telescope dRobberto et alj2010h . 
fully covering the area imaged with the HST. In particular, for 
our analysis we selected the WFI data, which were taken in the 
U, B, V and I passbands, in a narrow-band filter centered on the 
H a line and in a medium resolution filter centered at 6200 A. 
This provides a R-band line-free equivalent passband for stars 
K-type and earlier, whereas for stars M-type and later it is cen- 
tered on a temperature sensitive TiO absorption band which al- 
lows a reddening ind ependent measurement of T e ff for cool stars 
dDa Rio et alJl2 009). WFI measures were used only for stars not 
present in the ACS dataset. To avoid introducing color errors due 
to source variability, we have not combined the ACS and WFI 
data, as they were obtained at different epochs. 

In order to reliably build the HR diagram one has to com- 
bine the photometric observations with spectral types. We have 
therefore down selected from our original sam ple of 2 7 93 ob - 
jects 807 stars with spectral type given bv lHillenbrandl (fl997b . 
100 stars with spec tral type f rom a follow-up to th e study of 
IStassun et all d 19991) (see al so iDa Rio et all d2010ab l and 150 
stars with spectral type from lDa Rio et alj d2009). obtained with 
narrow -band photometry at 6200 A on a TiO absorption feature. 
In total we have therefore 1057 stars with known spectral type 
and photometry in at least three bands at optical wavelengths. 
The completeness of our spectroscopic sample will be analyzed 
in Sect. EH 



3. HR diagram 

3.1. Construction 

The construction of an accurate HR diagram requires the knowl- 
edge of several quantities and transformation relations. First, the 
spectral type has to be converted to effective temperature. For 
this purpose we used the rel ation between s pectral type and ef- 
fective temperature given bv lLuhmarj dl999l) . which is tuned to 
bring the presumably coeval components of the young system 
GG Tau on the same isochrone, according to the iBaraffe et al.l 



dl998h models. For M-type PMS stars this relation is interme- 
diate between dwar fs and giants, and as it has been shown by 
IDa Rio et al.l d2009 ), it represen t s an i mprovement over the re- 
lation assumed by Hillenbrand] dl997l) . Using this transforma- 
tion the temperatures of our sample range from ~ 2500 K to 
~ 8000 K. Second, to estimate the bolometric luminosity of 
each source one has to convert the observed magnitudes to ab- 
solute magnitudes adopting values for extinction, distance mod- 
ulus, and bolometric corrections. One should also account for 
non-photospheric sources of luminosity like, e.g., accretion lu- 
minosity. 

For the ACS data we used as a reference the F850LP filter 
(z-band equivalent) using the following relation: 

log(L*/L G ) = 0.4- 
[zacs,g - (zacs,* - Azacc - A z + BC Z ^(T) - DM)} 

where Zacs.o is the absolute magnitude of the Sun, Zacs,* is the 
magnitude of the target star, Az acc is the estimated contribution 
to the observed brightness due to disk accretion, A z is the ex- 
tinction, SC.- * is the bolometric correction and DM the distance 
modulus. For the WFI data we used an analogous at z expression 
for the / band. 

To estimate the various terms we computed synthetic pho- 
tometry on families of stellar atmospheric models using the 
SED fitting package Chorizos dMaiz-Apellanizll2004l) . The x 2 ~ 
minimization SED-fitting enables us to asses the correspondence 
between observed and synthetic photometry up to N=4 physical 
parameters. In other words, Chorizos solves the problem of find- 
ing which model SED is the most compatible with the observed 
colors in a model N-dimensional parameter space. We as sem- 
bled spectral libraries from various authors: iKuruczl rfl99l for 
effective t emperatures larger tha n 8000 K, NextGen model atmo- 
spheres dHauschildt et al.lll999b for 500 K < T < 8000 K and 
AMES-MT 2000 models dAllard et alJl2000h for T < 5000, all 
of th em with solar metal licity, following Da Rio et al.l d20 10a). 

IDa Rio etail d2010al) show that colors of ONC members dif- 
fer from those of main-sequence dwarfs in the M-type regime as 
well as the atmospheric models do not accurately match the pho- 
tometry of the ONC . Therefore, we limited the parameter space 
to models correspondin g to the empirical calibration suggested 
bv lDa Rio et al.l (I2010al) as a good approximation for the ONC. 
This empirical isochrone was built to obtain a good match for a 
selected sample of non-accreting low-extinction members of the 
ONC, using the ZAMS intrinsic colors for T > 3700 K, the 2 
Myr synthetic colors for T < 2900 K, and a linear interpolation 
between the two for 2900 >T< 3700 K. 

Moreover to account for the accretion luminosity that 
may affect our young sources, we added a typical accre- 
tion spectrum to the fami ly of template spectra. According to 
ICalvet & Gullbringl J 1998b . the total accretion column emission 
can be well approximated by the superposition of an optically 
thick emission from the heated photosphere below the shock 
with an optically thin emission generated in the infalling flow, 
in a 3:1 ratio. The optically thick component can be modeled 
by a black body at T e ff = 6000 - 8000 K, and we assumed 
T e ff — 7000 K. For the optically thin emission we used Cloudy 
dFerland et al.ll 998) to generate the spectrum of an optically thin 
slab with n = 10 8 cm 3 also at T =s 7000 K . We added the 
combined spectrum to our reference stellar models with ratios 
of \og{L accr I L tot ) increasing from -3 to 0. 

We ran Chorizos for each source of our sample, assuming the 
metallicity and the stellar effective temperature as fixed parame- 
ters and two different reddening laws, corresponding to the stan- 
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Fig. 1. Extinction and Accretion Luminosity. Top: Distribution of the 
extinction for the sample of 1057 stars from the ACS and WFI cata- 
logs for which we have estimated extinction and accretion luminosity. 
Bottom: Distribution of the accretion luminosity for the same sample of 
objects. The accretion luminosity is expressed in terms of L accr /L m and 
in logarithmic scale ranges from -3 to -0. 



dard Ry = 3.1 and Ry = 5. 5, by some considered more appro- 
priate for the Orion Nebula ( John son! 1967tlc"ostero & Pei mbert 
1970). Therefore, extinction and accretion represent the vari- 
ables of our fitting procedure. Chorizos returned the most likely 
value of extinction and accretion luminosity for each source, for 
each reddening law. With Chorizos we can fit all available pho- 
tometric data at the same time to deri ve extinction and acc retion 
luminosity for each source, whereas iDa Rio et alJ (1201 Oal) con- 
centrated on 2 colors. Chorizos is therefore more robust against 
the occasional spurious data points, or against the data with the 
highest photometric error. In Fig. Q] we show the distributions 
of the extinction and accretion luminosity for the two choices of 
reddening law. For the standard reddening law, that we will adopt 
as main reference, the extinction distribution has a median value 
at Ay = 2.2 mag and shows that roughly 70% of the sources 
are extincted by less than 3 mag. Concerning the distribution 
of (Laccr/Ltot), we find an average value for the relative fraction 
of 5% and more than half of the stars in the sample showing 
an excess attributed to accretion w ith L accr /L tot greater t han 3% 
(comparable to what was found by [Da Rio et al.l (1201031) 1. 



The bolometric corrections appropriate for our population 
were derived directly from the integration of the synthetic spec- 
tra in z as well as across the whole range of wavelengths: 



BC(T) = -2.5 log 



J A S A F A (T)dA J A S A F A (Q)dA 



■/■- 



JT F A (T)dA JT F A (Q)dA 



(1) 



where F A (T) is the spectral energy distribution of the star of tem- 
perature T, while F A (&) is the solar spectral profile, and S A is the 
passband of the filter multiplied by the quantum efficiency func- 
tion of the detector. According to EquationQ] the bolometric cor- 
rection of the Sun is zero in each band, whereas the zero-points 
in this case are irrelevant as they cancel out in the calculation of 
the luminosity. 

We also integrated the fluxes of a NextGen synthetic solar 
spectru m, assuming the canonical value of My^ (V Johnson) = 
+4.83 dBinnev & Merrifield 1998). The absolute magnitudes of 
the Sun in the z and / bands of ACS and WFI turn out, respec- 
tively, M z = 4.002 and M/ = 4.020, computed with synthetic 
photometry on solar spectrum. Concernin g the distance value , 
we assumed d = 414 + 7 pc according to Menten et al] d2007l) . 
This corresponds to a distance modulus DM = 8.085. 

3.2. Results 

Having estimated the temperature and luminosity for each 
source, we produced the HR diagrams for the two values 
of reddening laws mentioned above. We then compared the 
position of the stars with three sets of evolutionary tr acks 
and isochrones: ISiess et al.l d2000h . iPalla & Stahlerl (1 19991) and 
iD'Antona & Mazzitellil dl997t) (extended in 1998). 

In Fig. [2] we show the HR diagram relativ e to the red- 
dening law Ry — 3.1 and ISiess et all (2000) calculations. 
Currently there is no family of evolutionary models showing a 
significantly higher consisten cy with dynamical mas s estimates 
(Hillenbr and & White! 120041) . Therefore we adopt ISiess et all 
(2000) calculations as our reference, as they provide the smallest 
apparent age spread on the HR diagram. Moreover, according to 
IDa Rio et al.l (l2010al) . these models reduce the systematic corre- 
lation of average age with mass, more prominent when using 
IPalla & Stahlerl d 19991) isoch rones. or ID'Antona & Mazzitellil 
( 1994 models dHillenbrandlll997t) . In the diagram the objects 



are indicated according to dataset and available photometry, to- 
gether with evolutionary tracks and isochrones. In the next sec- 
tion we will analyze the accuracy of our estimate of the luminos- 
ity to quantify how much of the observed scatter can be due to 
observational uncertainties. Then, in Sect. |5j we will explicitly 
derive the age of each star through interpolation over the tracks 
and the isochrones to study the age distribution of our sample 
and to quantify the amount of scatter in terms of age spread. 

4. Sources of uncertainties 

As mentioned before, there are concerns regarding the reliabil- 
ity of derived ages, as well as masses, from the HR diagram of 
young clusters, due both to the discrepancies between different 
PMS evolutionary models and to the various uncertainties af- 
fecting observations. Whereas finding the optimal model goes 
beyond the scope of this wo rk, we can take a fresh look at the 
observational errors listed by Hart mannl d2001l) and discuss how 
they may affect our results on the ONC. First of all we present 
the analytical estimate of the uncertainties on luminosity and 
age, whereas in Sect. 14.21 we describe a Monte Carlo approach 
for the estimate of the same quantities. 
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Fig. 2. HR diagrams for the Orion Nebula Cluster. Black symbols represent ACS sources: circles are the stars with four-bands photometry while 
squares represent stars with complete three-bands photometry. Red squares indicate respectively WFI objec ts with 4 bands photometry, while red 
diamonds the WFI stars with three-bands photometry. Superimposed are the evolutionary calculations bv lSiess et al.l feOOO), Solid lines represent 
isochrones, the thickest one indicate the 1 Myr isochrone. Dotted lines represent evolutionary tracks and they range from 0.1 up to 5 M Q . 



4.1. Analytical estimates 

Distance. The VLBI (Very Long Baseline Array) distance ob- 
tained by measuring the trigonometric parallax of radio sources 
in the Orion Nebula Cluster gives d = 414 + 7 pc. The stan- 
dard deviation is larger than the transversal size (about 3 pc) 
of the cluster, and therefore still dominates over the intrinsic 
spread of distances within the cluster. We can therefore assume 
cr(logL) = 0.015. 



Sp ectral type. The spectral types given by iHillenbrandl 
(1997) are accu rate to within one sub-type. According to 
iHartmannl (120011) errors on the spectral type tend to cancel; e.g a 
star which is erroneously assigned an earlier spectral type would 
have a larger age at a constant luminosity, but the extinction 
would be also overestimated with consequent erroneously in- 
creased luminosity . Aver aging over the spectral type range of 
K6-M 1 . IHartmannl d200 ll) found that an error within a sub-type 
in spectral type corre sponds to an error on age of cr(log t) = 0.02. 
iDaRioet al.1 (l2010al) has shown that even when varying the ef- 
fective temperature within its error and consequently changing 
the bolometric correction in order to calculate the luminosity, 
due to the variation in extinction, the error on the luminosity is 



still par allel to the isochrone (see e.g. Figure 17 in lDa Rio et al.l 
d2010il) ). 

Extinction. Chorizos provides a typical uncertainty on the 
extinction (cr(Ay) ~ 0.2 mag). Using equation Q] to propagate 
errors, and knowing that A/ ~0.5 Ay and A z ~0.4 Ay, this trans- 
lates into cr(logL) -0.05. 

Photometric Variability. According to iHerbst et al.1 (|2002) 
the average variation in magnitude due to photometric variability 
for the ONC is Al = 0.18 mag. Following the same authors, we 
use the standard deviation to characterize the level of variability, 
since a general peak-to-peak measure can be affected by sin- 
gle accidental error. This standard deviation (0.075 mag) gives a 
logarithmic uncertainty on luminosity of <x(log L) = 0.03. 

Unresolved Binary companions. With a pixel scale of 50 mil- 
liarcsec, ACS is able to resolve sources of comparable brightness 
with separation larger than ~ 50 AU at the distance of the clus- 
ter. Assuming the standard lognormal peri od/separation distribu- 
tion t hat peaks at approximately 30 AU (Duquennov & Mavor 
1 19911), and a binary fractio n similar to the one observed in the 
field (Re ipurth et al.l 120071) . we have a probability of 20% of 
unresolved binary com panions. With th is probability the ba- 
sic model described in IHartmannl (1200 ll) allows us to derive a 
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logarithmic error on luminosity due to unresolved binaries of 
cr(logL) = 0.016. 

Accretion. Our analysis allows us to derive the amount of 
accretion luminosity that affects each source and to subtract 
it. However, the error on L acc is given by Chorizos in terms 
of cr(\og{L accr I L tot )), and consequently depends on <x(logL rof ) 
which is what we want to determine. The best way to estimate 
the effect of this uncertainty on the luminosity is through Monte 
Carlo simulations. Therefore, we neglect this contribution in this 
preliminary calculation. 

The quadratic sum of all these errors gives an overall un- 
certainty cr(logL) = 0.06, w hich is roughly 1 /3 of the value 
cr(logL) = 0.16 reported by iHartmannl (1200 ll) for the Taurus 
cloud. This reflects the better quality of our data set. 

It is known that the ages of low-mass pre-main-sequence 
stars on the Hayashi tracks exhibit an approximate power-law 
relation t o the stellar lum inosity L for each mass, of the form 
t oc Zr 3 / 2 dHartmannl2009l) . Because of this relationship between 
luminosity and age, the error distribution in age takes the form 
5t/t ~ 6\ogt, with 



cr(log*) = (3/2)cr(logL). 



(2) 



Therefore, from the uncertainty on the luminosity one can de- 
rive the error on the age, which turns out to be <x(logf) = 0.09. 
As shown also by Hartmannj (1200 ll) . errors of this order may 
preserve a time dependence of the star formation rate possibly 
present in the observations, a dependence which instead is lost 
when the observational errors are cr(log f) >0.2. 

Since our preliminary estimate of the overall uncertainty 
turns out to be smaller than the limit indicated by Hartmann, we 
can be fairly confident that, from the observational point of view, 
our uncertainties may allow us to infer real age trends from the 
analysis of the HR diagram. However, to understand to which 
extent this can be done, a more complete analysis is needed (see 
Sect.©. 

4.2. Monte Carlo simulations. 

After having determined a rough estimate of the uncertainty on 
luminosity, to better quantify the error on our age estimates, we 
have run Monte Carlo simulations to evaluate the dispersion in 
the estimated absolute luminosity of each object due to all the 
sources of uncertainty illustrated in the previous section. By gen- 
erating artificial errors for each star we can more accurately de- 
termine the overall uncertainties on the the bolometric luminos- 
ity. 

To evaluate the extinction and the accretion luminosity we 
have run Chorizos varying the photometry according to the pho- 
tometric errors and adding variability , following the var iabil- 
ity distribution in the I band given bv lHerbst et all (l2002h . We 
computed the amount of variation in the other bands through 
synthetic photometry on artificial stars with standa rd spot tem- 
peratu res, radii and filling factors for T Tauri stars dBerdv ugina 
2005). Each time we have also varied the value of the distance 
within its error and added a shift in luminosity due to binarity, as- 
suming the probability of unresolved companions given in Sect. 

s 

We ran N=200 trials for each source, determining each time 
the error in the estimate of the luminosity and averaged the dis- 
persion over the whole sample, obtaining values for <x(logL). 
We assume the value of cr(logL) = 0.10 as representative of the 
mean uncertainty we have on the luminosity. 

According again to Equation|2] this error translates in an un- 
certainty on the estimate of the ages of cr(logf) = 0.15, which 




Siess (2000) 
Mean - 6.36 
Sigma - 0.38 



_ Rv - 5.5 



Mean = 6.30 
Sigma = 0.46 




log(age/yr) 

Fig. 3. Distributions of stellar ages. Histo gram of the num ber of stars 
respect to the stellar ages, as derived from lSiess et alJ f2000) evolution- 
ary calculations for the reddening laws R v = 3.1 (top) and R v = 5.5 
(bottom). 



is larger than what we derived from basic error propagation 
(cr(log f) = 0.09). A larger error is expected, since we are con- 
sidering also the uncertainty in the estimate of the accretion lu- 
minosity, that was not kept into account in the previous calcula- 
tio n. Still, thi s value is below the cr(log t) = 0.2 threshold given 
bv lHartmannl yOOl), and clearly remains much smaller than the 
observed scatter (see Sect.O. 



5. The age distribution 

In Fig. [3]we show the logarithmic distribution of stell ar ages ob- 
tained for the two reddening laws we adopted and the Sies s et alj 
(2000) evolutionary models. The average ages and age disper- 
sions for the ONC that we derived with this set of evolutionary 
tracks and for the different reddening laws are listed in Table [2] 

The cluster age turns out to depend on the reddening law. 
In general, ages derived assuming a reddening law of Ry = 3.1 
are larger than those derived assuming Ry = 5.5. This is ex- 
pected, since high values of Ry imply higher extinction values 
and therefore higher intrinsic stella r luminosi ty and younger age. 
The mean age predicted bv lSiess et alj d2000l) for the sample is 
~ 2.3 Myr and ~ 2 Myr for Ry — 3.1 and Ry = 5.5, respectively. 

Concern i ng va riability, using the information supplied by 
iHerbst et al.l (120021) . it is possible to select a sub-sample of 508 
objects showing a peak-to-peak variability AI < 0.4 and 298 
objects with AI < 0.1. For these two samples the total error in 
the logarithm of the luminosity, as estimated through basic error 
propagation, goes to cr(logL) = 0.053 and cr(logL) = 0.051, 
respectively (for the complete sample cr(log L) = 0.06, see Sect. 
[4}. In Fig. 3] we show the distributions of the ages for the differ- 
ent samples and in Table [3] we summarize the results. A direct 
comparison with the estimates obtained for the full sample does 
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Table 2. Mean age and age dispersion 





Rv = 


3.1 


Rv = 


5.5 


PMS model 


< log(age) > 


a log(age) 


< log(age) > 


crlog(age) 


Siess et al. (2000) 


6.36 


0.38 


6.30 


0.46 


D'Antona & Mazzitelli (1998) 


5.93 


0.59 


5.82 


0.71 


Palla&Stahler(1999) 


6.14 


0.49 


6.07 


0.53 




Fig. 4. Distributions of st ellar ages for low variability stars. The stel- 
lar ages are derived from Siess et al. ( 2000) evolutionary calculations 
for the reddening laws R v = 3.1 (top) and R v = 5.5 (bottom). The 
open histogram represents the total sample. The yellow area indicates 
the sam ple of 608 stars for which we have information about variabil- 
ity from Herbst et al. (2002). The red hatched histogram represents the 
sources selected with Al < 0.4 mag while the black one the stars with 
A/ < 0. 1 . The average age of the cluster and the age spread are relative 
to the third sample. 
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Fig. 5. HR diagrams for the Orion Nebula Cluster. Same legend as in 
Fig. [2] From the top to the bottom are superimpose d evo lutionary tracks 
and is ochrone from lD'Antona & Mazzitelli dl997h and lPalla & Stahleil 
( 1999). Dotted lines represent evolutionary tracks and they range from 
0.017 to 3 M Q for D'Antona & Mazzitelli and from 0. 1 to 6 M e for Palla 
& Stahler. 



not show substantial differences both in the average age and age 
spread. The Kolmogorov-Smirnov test (KS test), which returns 
the probability that two distributions were drawn from the same 
parent sample, confirms that the samples are consistent with the 
original population, ie. variability does not affect significantly 
the mean and the spread of the age distribution. 

As mentioned before, we have considered three differ- 
ent families of PMS e volutio nar y tracks a nd isoch rones: 
D'Antona & Mazzitellil dl997t) . ISiess et all (120001) and 
Palla & Stahlerl (1 19991) . In Figs S and we present the HR 
diagram and the age distribution according to these different 
sets and in Table|2]the mean ages and the age spreads predicted 
for the Orion Nebula Cluster. The differences in the estimate 
are clearly visible. First of all, this discrepancy affects both the 
mean age of the cluster and the age of th e single star. Regardless 
of the value of the reddening parameter. iD'Antona & Mazzitellil 
(1997) consistently predicts younger ages, which appear to be 
consistent with the results obtained for the same evolutionary 
calculations by Hillenb rand! (1 19971) . Second, the value of the 



mean age depends also on the predictive capability of the model 
for this particular cluster. Since the different sets cover slightly 
different areas in the HR diagram (Fig. [5}, certain stars may or 
may not be left out of certain models. Finally, a relative age 
spread appears to be robust regardless of which family of models 
is used, although the am ount of this s pread is model dependent. 
The scatter predicted by ISiess et al . (2000) is the smallest one 
with both the standard and the anomalous reddening law. 

In conclusion, whereas from the observational point of view 
our sample enables us to place the stars in the HR diagram with 
fairly high precision, different PMS evolutionary models provide 
significantly different ages. But regardless of the choice of the 
set of models, there appears to be an age scatter which cannot 
be accounted for by observational errors. Star formation in the 
ONC has mainly occurred within the last 10 Myr, with ~ 70% 
of the stars younger than 3 Myr. There is no evidence of a single 
star formation burst as the distribution spans a wide range of 
ages with a scatter of a few Myr s, in agreement with what was 
found bv lHillenbrandl(ll997l) and ba Rio et al.l d2010al) . 
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Herbst Sample 


First selection 


Second selection 


Rv 


<log(age)> crlog(age) 


< log(age) > crlog(age) 


<log(flge)> crlog(flge) 


3.1 


6.33 0.36 


6.31 0.34 


6.31 0.36 


5.5 


6.37 0.43 


6.25 0.42 


6.27 0.41 




Fig. 6. Distributions of stellar ages for different PMS evolutionary 
models. These histograms are relative to the reddening laws R v = 3.1 
(top) and R v = 5.5 (bottom). The y ellow histogram represents the 
results obtained with interpolation on iD'Antona & Mazzitellil l | 19971) 
isochrones. The red hatched histo gram is relative to Siess et al. (2000) 
while the black hatched one shows Palla & Stahlej i 19991) predictions. 



5. 1 . Completeness Analysis 

The HRD and the age distribution presented here, are relative to 
the subsample of ONC objects with at least three-band photom- 
etry and known spectral type. This represents only a fraction of 
our entire photometric catalog of the ONC, and its completeness, 
which worsen towards fainter luminosities - or lower masses - 
is basically depending on the availability of spectral types, and 
not on the photometric detection. This is because, whereas our 
HST photometry extend s well into the BP m ass range and goes 5 
magnitudes fainter than Hille nbrand! d 1 9971) . spectral types reach 
at most the lowest masses in the stellar regime. 

Towards this limit, the subsample becomes also extinction 
limited, in the sense that very-low mass members will have 
an available T e s estimate only if their Ay is not significantly 
large. The reason for this comes from source of our spectral 
types: the majority of these are from the optical spectroscopy of 
Hillenbrand (1997); in that work optical spectra were obtained 
selecting randomly the candidates, but only for sources with pre- 
viously available V and / photometry, and clearly were limited 
to lower mag nitude limits. Other spectral types were derived by 
iDaRioet al.l f2009) using a photometric index computed from 
ground-based measures in the V and / bands, plus a medium- 



band filter centered on the TiO feature at A ~ 6200A; this second 
set reaches spectral types as late as M6, for moderate extinction. 

Besides the optical lu minosity constraint, n o ad ditional crite- 
rion w as imposed (both in Hillen brand! d 1997b and lDaRio et al.l 
(2009)) for the selection of sources with assigned spectral types. 
Consequently, the only selection effect which affects the com- 
pleteness of our subsample is the apparent V magnitude of the 
members, since due to the intrinsic color and reddening of typ- 
ical ONC populations the detection is generally easier at longer 
wavelengths. 

Therefore, we computed the completeness distribution of our 
subsample as a function of V magnitude. Even at the faintest 
end of range of V magnitudes, the ACS photometric catalog 
can be considered complete reaching several magnitudes deeper. 
Therefore, the ratio of the two V-band luminosity functions (in 
the subsample and for the whole ACS catalog) provides automat- 
ically the V-band completeness function. More precisely, this 
computation was done separately in three subcases: the ACS 
stars, WFI stars in the FOV of ACS, and the WFI stars outside 
the ACS FOV. Since at the bright end of the LF, the ACS cata- 
log is deficient of sources with respect to WFI, due to saturation, 
we complemented the ACS catalog with the missing bright WFI 
sources falling in the same area. By adding WFI sources outside 
the ACS FOV, we are allowing the completeness function to be 
> 100% for some bright V-magnitudes; however, this is not a 
problem, since any multiplicative re-normalization of the com- 
pleteness function is absolutely arbitrary and irrelevant for our 
purposes. 

We want now to pass from the V-band completeness to the 
completeness as a function of stellar parameters, like T e g and 
log L, or better, mass and age. This cannot be done directly, be- 
cause differential reddening, as we said, introduces an extinction 
limited bias. The computation is still possible in a probabilistic 
way, through a monte-carlo simulation. The method is identical 
to that used in lDa Rio et al.l d2010al) : we consider a pair of values 
of mass and age. Assuming the evolutionary models, we convert 
this point in T e s and log L. Then we apply backwards bolomet- 
ric corrections and intrinsic colors to covert this point further to 
the V vs V - / CMD, assuming Ay — 0. Having the intrinsic Vo 
magnitude of the (mass, age) point, we redden this test star ran- 
domly, with values of Ay drawn from the reddening distribution 
of the ONC, computed on the brightest half of the population. 
Each of this "reddened" test stars will have a different (fainter) 
V magnitude than Vo; and therefore a different completeness. By 
simply averaging all these completeness ratios, we derive the ex- 
act completeness of the point in the mass-age plane considered at 
first. Iterating this procedure for each (mass, age) pair, we derive 
a 2-dimensional completeness function in the mass-age plane. 

Figure [7] shows the completeness contours in the mass-age 
plane for the two different reddening laws, whereas in Figure [8] 
we pr esent the complete ness-corrected age distributions accord- 
ing to lSiess et al.l {2000) evolutionary models for the two differ- 
ent reddening laws. We observe that completeness decreases to- 
ward lower masses and higher ages, therefore the completeness 
correction does not change significantly the average ages and the 
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Fig. 7. Contour plots of the completeness function in the mass/age 
plane, com puted from our simulations for the evolutionary models of 
ISiessetalJOOOOh and for R v =3. 1 (top) and R v =5.5 (bottom). 
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Fig. 8. Completeness-corrected age distribution. The hatched his- 
togram represents the age distribution of our sample. The open his- 
togram shows the age distribution of the completeness-corrected sample 



measured age spread (< log(age) >=6.40 and cr log(age)=0.49 
for Ry=3.1) but increases the old age end of the distributions. 

However, optical photometry alone is often not enough to 
identify the complete stellar population in young clusters. Near- 
infrared surveys of the ONC in the past decades have revealed 
the presence of a large number of highly extincted sources in 
this r egion (e.g. iHillenbrand & Hartmannl[l998t iMuench et alJ 
12002V even though a more recent near-infrared study has shown 
that a majority of st ars (-70%) of the ONC has Ay <5 
( Rob berto et al.l d2010h . Scandariato et al. in prep.). Moreover 
the contamination becomes an important issue with increasing 
sensitivity of the surveys. Therefore it is not correct to assume 
that everything that has been observed in the near-infrared is 
a cluster member. In any case, because our survey is not com- 
plete down to 0.25 M for more deeply embedded associa- 
tion members, we have tested whether our sample is biased to- 
wards older stars. We considered a mass-age-Ay limited sample 
(M<0.25M G , t < 10 7 Myrs and A v <5) of 647 stars in total. We 
divided this subsample in two halves with equal number of ob- 
jects according to low (Ay <2) and high (Ay >2) extinction. We 
have used the KS test to check whether the age distributions of 
the two samples are drawn from the same parent population and 
we are unable to reject the hypothesis that the distribution are the 
same (P(d)fcs =0.7 6). Since we do not see evidence for a trend of 
age with increasing Ay, we argue that having most of the stars, 
though not all, in our sample allows us to study the age distribu- 
tion of the Cluster based on a lower extinction sub-sample. 



5.2. Membership 

It could be argued that some of the scatter observed in the HR 
diagram and the age distribution might be due to non-member 
contamination. The Orion complex itself consists of three sub- 
groups, the ONC (Ori on subgroup I d) an d other three subgroups 
of stars (la, lb, and Ic. lBrown et al.ldl994h V The most likely sub- 
group from which we would see contamination in our data is 
the Orion subgroup Ic, since it is located along the same line 
of sight of the ONC. Orion Ic is t hought to be less than 5 Myr 
old according to lBrown et al.l (1 19941) . thus we would expect little 
differences between the isochrones of the two subgroups. 

In any case, for the ONC objects also present in H97 we 
utilized the membership given in that work, which was take n 
from the proper motion information of I Jones & Walked (1 19881) . 
to study the age distribution of the confirmed members. Among 
our 1057 stars, 701 have proper motion measurements and are 
thought to be ONC members with a probability P >50%. In Fig. 
[9] we show the age distribution of the whole sample (open his- 
togram) and of the 703 co nfirmed members (hatched histogram) 
in the case of ,Ry=3.1 and lSiess et al.l (2000) evolutionary mod- 
els. The KS test give a probability of 17% that the two distribu- 
tion are drawn from the same parent population and the spread 
(cr log(a,ge)=0.40) obtained for the subsample does not differ 
significantly from cr log(a,ge)=0.43 derived for the whole sam- 
ple. 

For these reasons we are confident that contamination does 
not contribute very much to the scatter observed in the age dis- 
tribution of the ONC. 
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Fig. 9. Age distribution of confirmed members. The open histogram 
represents the age distribution of the complete sample. The hatched his- 
togram shows the age distribution for the stars also present in the H97 
catalog with me mbership probab i lity P >50% according to the proper 
motion survey of lJones & WalkeJ d 1988h . The KS test gives a probabil- 
ity of 17% that the two distributions are consistent. 



6. Age spread 

The observed dispersion in age represents an upper limit to the 
real age spread, being a combination of true age spread and scat- 
ter due to observational uncertainties. These, however, cannot 
explain the scatter of ~ 1 dex observed in the HR diagram, as 
we have found in the previous sections. To quantify, or at least 
to put a constraint, on the real age spread it is necessary to run 
a second Monte Carlo simulation accounting for observational 
errors. 

In order to increase the reliability of our results, we lim- 
ited the analysis in the mass range in which our completeness 
is above 80%. Our first step has been to test the hypothesis of 
an exactly coeval cluster, simulating a sample of stars with lumi- 
nosities and temperature corresponding to a 2.2 Myr isochrone 
( mean age va l ue for the ONC estimated in Sect. [3j according 
to ISiess et aD (l2000t) evolutionary model. We have then varied 
their luminosities assuming a Gaussian distribution of errors 
with cr(logL) = 0.1 (see Sect.H), rec overed their ages from in- 
terpolation over the Siess et al. (2000) isochrones and evolution- 
ary tracks and studied the overall age distribution. The spread 
around the isochrone is cr(\ogf) = 0.13. We compared this dis- 
tribution with the observed one using the KS test and we found 
a probability of 10~ 19 that the two distributions are drawn from 
the same parent population. The comparison is shown in Fig.fTUl 
This formally confirms that the observational uncertainties alone 
are not able to account for the age spread observed in the region. 

We then performed other simulations to find the amount of 
age spread needed to reproduce the observations, in addition to 
the dispersion due to observational uncertainties. In Fig. fTTIwe 
show the comparison with the simulations in the cases of spreads 
of 1, 1.5, 2, 2.5, 3, 3.5 and 4 Myrs. The probabilities given by the 
KS test are also shown in each panel of Fig.QT] Since for spread 
larger than 3.5 and smaller than 1.5 Myr the KS test probability 
drops below 1%, we suggest the real age spread to be within 1.5 
and 3.5 Myr. 

Therefore, on the basis of the isochrones, we find a spread 
of roughly 2-3 Myr in age, which is in good agreement with 
the results presented by iDa Rio et al. (2010a), and other tech- 
niques for inferring stellar ages also confirm the same results 
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Fig. 10. Test of coeval population. (Top) The open histogram represents 
the age distribution simulated f or a sample stars uniformly distributed 
along a 2.2 Myr isochrone from lSiess et 23(2000), assuming an obser- 
vational uncertainty on the luminosity of cr(log L) = 0. 10. The spread in 
age is <x(logf) = 0.13. The hatched histogram shows the observed age 
distribution for the observed completeness-limited subsample. (Bottom) 
Cumulative functions of the the two distributions, solid and dotted line 
represent the observed and simulated sample respectively. 



dMuench et al.ll2008h . Another way to interpret this result is that 
we observe a scatter in stellar radii at any given mass. The same 
evidence was found bv lJeffriesI (|2007) using the rotation periods 
and project ed equatorial ve locities (v sin z) for a sample of stars 
in the ONC. Jeffries] d2007l) studied the distribution of R sin i and 
modeled it as dispersion of true radii, inferring a spread in stellar 
ages. This spread is again approximately a few Myrs. 

Wit h refer ence to the recent at t empts , iBaraffe et all j2009), 
iJeffriesI (120111) and iLittlefair et aD (1201 ll) with different meth- 
ods suggest that the spread in radius observed in different star 
forming regions may be due to the accret ion history and not to 
a real spread in age. IBaraffe et al.l d2009l) present a theoretical 
calculation for t he expected sc atter in luminosity due to episodic 
accretion, while IJeffriesI (1201 ll) find that stars with different ages 
do not show differences in accretio n and disk excess an d propose 
they have the same age. Finally. ILittlefair et al.l (1201 ll) observe 
a correlation between rotation velocity and luminosity, which is 
opposite in respect to what is expected in the case of a real age 
spread. This, according t o the authors, might be supportive of 
the sc enario presented in IBaraffe et al.l d2009l) . iHosokawa et al.l 
(1201 ll) . on the other hand, show that the broadening of the PMS 
sequence in the HR D due to episodic accretion, as computed by 
Bara ffe et al.1 (120091) . was significantly overestimated, indicating 
that even this scenario might be incomplete. 
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Fig. 11. Age spread test. (Top, from left to right) The open histograms represent the age distribution simulated for a sample stars with mean age of 
2.2 Myr, assuming an observational uncertainty on the luminosity of cr(logL) = 0.10 and with intrinsic age spreads of 1, 1.5, 2, 1.5, 3, 1.5 and 4 
Myr respectively. The KS test probability between the two distributions is also given. (Bottom) Cumulative functions for each simulation. We use 
the same legend as in FigllOl 



At the moment we do not exclude a superposition of a real 
age spread and episodic accretion, especially since a model 
quantifying whether the accretion-induced luminosity spread is 
able alone to account for the observations has not been devel- 
oped yet. In any case a more complete scenario still needs to be 
investigated. 

We remark, finally, that the Orion Nebula hosts a few sites 
of ongoing star formation, like e.g. the BN/KL system, Ori-1 
South, and an extremely dense filament of gas and dust stretch- 
ing through the cen ter of the ONC whic h hosts at least 45 proto- 
stellar candidates (Mue nch et al.l l2008'). In our study of the op- 
tically visible population we could not taken into account these 
embedded regions. Our finding that the mean age of the ONC 
is 2 Myr and that the star formation appears to have occurred 
between ~l-3 Myr ago is consistent with the currently detected 
level of ongoing star formation. However, if future observations 
reveal dozens more deeply embedded protostars over the same 
field of view as we observe several hundreds of T Tauri stars 
with ages of 1-3 Myr, the overall picture might change. 

7. Spatial age distribution 

In Sect. [5]we derived the age distribution for the ONC according 
to different evolutionary models and reddening laws. We discuss 
now whether this distribution holds true at different location in 
the cluster. 

We limit again our investigation of the cluster age structure 
to the subsample of stars for which our completeness is above 
80%. In Figure [12] we show the spatial distribution on the sky 
of this subsample of ONC stars, binned by age. Each bin has 
roughly the same number of stars. As can be easily seen, the 
youngest stars are more clustered towards the center, while the 
oldest ones are almost homogeneously distributed i n space. This 
eviden ce for age segregation was already noticed bv lHillenbrandl 
(fl997h . A two dimensional KS test also confirms this first-look 
impression. According to the KS test, whereas comparing the 
spatial distributions of stars younger than the mean age of the 
cluster (first two panels) we get a probability of 16% that they 
are consistent with each other, we find a probability <1% that 
they are in agreement with the spatial distributions of the stars in 
the last two bins. Moreover, more than 65% of the stars younger 
than 1 Myr are located within one parsec from the central O star, 
l Ori C, commonly assumed as center of the cluster. 

This preliminary result suggests a possible trend of ages with 
spatial distribution. However, a part of this trend might be biased 
by spatially variable incompleteness. In fact, while our com- 



pleteness function is representative, on average, of the whole 
stellar sample we considered, the completeness varies with po- 
sition. Specifically, it is worse towards the center of the system 
where the sky background is brighter. In any case, if this "age 
segregation" is confirmed, this may have implications for the 
study of the dynamical evolution of the cluster. 

8. Dynamical evolution of the ONC 

In the past years there have been man y attempts to understand 
the dy namical histor y of t he ONC (e.g. lHillenbrand & Hartmannl 
1998; lAllison et al.ll200^). In particular the evidence of mass 
segregation ( Hil lenbrand! ! 1997; Hill enbrand & Hartmannl 1998) 
in such a young cluster highlights the problem of distinguishing 
between primordial and dynamical mass segregation. In fact, the 
cluster crossing time is ~ 2 Myr, roughly the same mean age of 
the cluste r, and an even longer time is r equired for the cluster 
to relax (Hille nbrand & Hartmannl[l998l) . Thus, mass segrega- 
tion could be primordial, meaning that the ONC formed with the 
most massive stars concentrated near the center. On the other 
hand, simulations of initially cool and clumpy clusters have 
shown that with these initial conditions mas s segregation coul d 
also be dynamical on very short timescales (Alli son et al.l l2009). 
However if the observed spread in radii can be interpreted as a 
true age spread and the star formation in the ONC lasted for at 
least 2 Myr, the dynamics for the youngest and oldest stars could 
turn out to be very different: the evolution of the oldest, earliest 
stars was dominated by the molecular cloud in which they were 
embedded, whereas the youngest stars in the clusters were sub- 
ject to the dynamical influence of the stars already formed. This 
highlights that models which attempt to explain dynamical evo- 
lution in this region are probably still incomplete. 

9. Conclusions 

In this work, taking advantage of the broad-band photometry ex- 
tracted from the large data-set of ACS/WFC observations col- 
lected for the HST Treasury Program on the ONC, we have 
derived the HR diagram for the Cluster, with the main goal of 
studying the age distribution and age spread in the region. 

Our method is based on the estimate of the extinction 
and accretion luminosity towards each source by perform- 
ing synthetic p hotometry on an em pirical calibration of atmo- 
spheric models (|Da Rio et al.l2010al) using the package Chorizos 
(Mafz-Apellanizi l2004i) . This tool allows us to derive accurate 
estimate of the intrinsic luminosity for each star. Taking into 
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Fig. 12. Spatial age distribution. The four images show the spatial distribution of stars within various age ranges. The youngest stars are clearly 
more centered towards the center than the oldest ones. 



account uncertainties on the distance, spectral type, extinction, 
unresolved binaries, accretion and photometric variability we 
can determine luminosity and age with logarithmic errors of 
cr(logf) - 0,10 a nd cr(logf) = 0.15, respectively. Compare to 
Hartmann (2001), we obtain errors small enough to preserve a 
time dependence of the star formation rate possibly present in 
the observations, which is instead lost when the observational 
errors are cr(log f) >0.2. 

The HR diagra m we derived show s a large spread in lumi- 
nosity. According to lSiess et al.l (12000) evolutionary calculations 
the mean age of the cluster is 2.2 Myr, and it is slightly depen- 
dent on the assumed reddening law. Much stronger is the de- 
pendence on the evolutionary models since choice of different 
isochrones and tracks affects both the mean age of the Cluster 
and the estimate of the age of the single star. Currently there are 
several evolutionary models available and our capability of in- 
dicating the best one is limited to secondary clues (e.g. spacial 
distribution and kinematics vs. average cluster age) and our prej- 
udice (e.g. lack of mass-age relation). However, the presence of 
an observed spread larger than the errors on the age is common 
to all the models. 



A preliminary study of the spatial distribution of the ages 
suggests age segregation in the ONC. However we cannot ex- 
clude the possibility of an apparent trend, due to a spatial de- 
pendence of completeness. A further investigation is certainly 
needed. 

To quantify the amount of intrinsic luminosity spread in 
terms of age spread, we ran also MC simulations taking into 
account the sources of errors mentioned before. First of all, 
the observations are not consistent with a coeval stellar popu- 
lation, even if we account for the apparent age spread due to 
observational uncertainties. Second, the observed age distribu- 
tion in the ONC is consistent with simulations of stellar popula- 
tion with intrinsic age spread between 1.5 and 3.5 Myrs, but not 
larger. These results confirm the presence of a real luminosity 
spread which cannot be accounted for only with observational 
uncertainties but reflects a real spread in the stellar radii for a 
given mass. We interpret this a true age spread in Orion, im- 
plying that star formation in this region lasted a couple of Myr. 
Some authors suggest that this spread might be due to the accre- 
tion history but at the moment the truthfulness of this scenario 
has not been proved yet. However any possible explanation for 
the spread observed in the HR diagram should reproduce this 
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amount of age spread. Within the current limitations, we con- 
clude that the age is 2 Myr and the age spread is ~ 2 Myr. 
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